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Abstract: On directed and undirected Barabasi-Albert networks the Ising model 
with spin S = 1/2 in the presence of a kind of noise is now studied through Monte 
Carlo simulations. The noise spectrum P(n) follows a power law, where P(n) is 
the probability of flipping randomly select n spins at each time step. The noise 
spectrum P(n) is introduced to mimic the self-organized criticality as a model 
influence of a complex environment. In this model, different from the square lattice, 
the order-disorder phase transition of the order parameter not is observed. For 
directed Barabasi-Albert networks the magnetisation tends to zero exponentially 
and undirected Barabasi-Albert networks , its remain constant. 
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Introduction 

This paper deals with Ising spin on directed and undirected Barabasi- 
Albert (BA) networks in the presence of a noise. Sumour and Shabat [H [2] 
investigated Ising models with spin S = 1/2 on directed BA networks j3] with 
the usual Glauber dynamics. No spontaneous magnetisation was found, in 
contrast to the case of undirected BA networks [H [5j [6] where a sponta- 
neous magnetisation was found below a critical temperature which increases 
logarithmically with system size. Lima and Stauffer [8] simulated directed 
square, cubic and hypercubic lattices in two to five dimensions with heat 
bath dynamics in order to separate the network effects of directedness. They 
also compared different spin flip algorithms, including cluster flips [9], for 
Ising-BA networks. They found a freezing-in of the magnetisation similar to 
[H [2] , following an Arrhenius law at least in low dimensions. This lack of 
a spontaneous magnetisation (in the usual sense) is consistent with the fact 
that if on a directed lattice a spin Sj influences spin Si, then spin Si in turn 
does not influence Sj, and there may be no well-defined total energy. Thus, 
they show that for the same scale-free networks, different algorithms give 
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different results. Recently Stauffer and Kulakowski [TU] simulated the Ising 
two-dimensional ferromagnet in the presence of a special kind of noise dif- 
ferent from the traditional temperature and Boltzmann probabilities, where 
the noise consist of at each iteration of a L x L square lattice with four 
neighbours for each sites, besides the above majority rule they select n times 
randomly a spin and flip it. The probability distribution function P(n) of 
these numbers n is taken as a power law, 

P(n) oc l/n a . (1) 
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Figure 1: Summed magnetisations versus time, T = 1, a = 1, and m = 2 for 
directed (BA) networks . 



In order to get this distribution, they determined random numbers r, 
homogeneously distributed between zero and one, and then took 

n = T I?r l ^- a \ (2) 
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Figure 2: Summed magnetisations versus time for 0.01 < T < 1.0, a = 1 
and m = 2 (left) and 0.1 < T < 1.0, a = 1 and m = 7 (right) for directed 
(BA) networks. 

for a < 1, where 1 < n < L 2 and 

n = Texp(rln(L 2 )) (3) 

for a = 1 and T determines the amplitude of the noise. Here we simu- 
lated the same model described above on directed and undirected (BA) net- 
works and our results are different from the results obtained by Stauffer and 
Kulakowski[10j. 

Results 

In Fig. 1 we show the dependence of the summed magnetisation M = 
versus time for a = 1, T = 0.01(eq.3) and connection number m = 2 
on directed (BA) networks with 500000 sites, when we started with all spins 
up. M relaxes exponentially towards zero. In Fig. 2 we plot same picture 
of Fig.l, but for several values of T and connection numbers m = 2 (left) 
and m = 7(right), where for smaller T the magnetisation tries to remain 
positive. In the Fig. 3 we plot the dependence of the summed magnetisation 
versus time for a — 1, several values of T(left), and for T = 10(right) even 
up to 10000 Monte Carlo steps with m = 2 on undirected (BA) networks 
with 10000 sites. The summed magnetisation remains positive and does 
not decay towards zero. It keeps fluctuating around a positive value of the 
magnetisation depending on the value of T. 
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Figure 3: Summed magnetisations versus time for 1 < T < 1000., a = 1 
and m = 2 (left) and for T = 10, a = 1 and m = 2 (right) for undirected 
(BA) networks. 

In Fig. 4 the magnetisation is shown as a function of time for one sample 
for T = 1, a = 1 and m = 2 on directed^ left) and un<izrecte(i(right)(BA) 
networks. For directed (BA) networks the summed magnetisation changes 
between positive and negative values with time, and for undirected (BA) net- 
works the magnetisation fluctuates around a positive value. In Fig. 5 we 
display the normalized magnetisation versus the temperature T for directed 
(BA) networks. It decays exponentially towards zero for 100 < t < 1000 
which shows that a phase transition for this model does not exist, in agree- 
ment with the results for Ising model on this network [TJ [2], but in disagree- 
ment with the results of Stauffer and Kulakowski[10] for square lattice. The 
absolute value of the magnetisation for some T values ( no shown here) fluc- 
tuates around a positive value of the magnetisation depending on the value 
of T, in agreement with Fig. 4. 

Discussion 

We have presented the simulation of majority rule disturbed by power- 
law noise on directed and undirected (BA) networks. We study this model 
for case of a = 1, eq.3, and values different from 1 (not shown here), eq.2, in 
both cases the results are the same. For directed (BA) networks the magneti- 
sation presents no phase transition and decays exponentially with time and 
with temperature. For the undirected (BA) networks the magnetisation also 
presents no indication of a phase transition in the presence of the noise; it 
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Figure 4: Magnetization versus time for one sample, m = 2 for directed (left) 
and undirected(nght) (BA) networks. 

assumes positive values that depend on the values of the temperature. These 
results are differentfrom the results obtained by Stauffer and KutakowskifTU] 
for a square lattice. 
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Figure 5: Normalized magnetization, average over 1000 samples, 10 networks 
and over 100 < t < 1000, m = 2 for directed (BA) networks. 
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